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We investigate the robustness of some recent results obtained for homo- 
geneous and isotropic cosmological models with conformally coupled scalar 
fields. For this purpose, we investigate anisotropic homogeneous solutions of 
the models described by the action 



with general F(<p) and V(4>). We show that such a class of models leads gener- 
ically to geometrical singularities if for some value of eft, F(4>) = 0> rendering 
previous cosmological results obtained for the conformal coupling case highly 
unstable. We show that stable models can be obtained for suitable choices of 
F((j)) and V((ft). Implications for other recent results are also discussed. 
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I. INTRODUCTION 



Recently P],[2| , we have investigated the dynamics of homogeneous and isotropic solutions 
of the cosmological model described by the action: 

S = J d 4 xy/^ {F(<p)R - d^d^ - 2V(<f>)} , (1) 

with F((p) = 1 — |0 2 , corresponding to the so-called conformal coupling, and V(<p) = 
Y<fi 2 — j 4 - Some novel dynamical behaviors were identified: superinflation regimes, a 
possible avoidance of big-bang singularities through classical birth of the universe from 
empty Minkowski space, spontaneous entry into and exit from inflation, and a cosmological 
history suitable for describing quintessence. Since one of the proposals of inflationary models 
is to describe our universe without finely-tuned parameters, these results would be viable 
only if they are robust against small perturbations in initial conditions and in the model 
itself. This is the question to be addressed here. 

We study the robustness of our previous results by taking two generalizations of the model 
considered previously: we relax the isotropy requirements (perturbations in the initial con- 
ditions) and we consider a general coupling F(<j)) (perturbations in the model parameters). 
Models with more general F(<j)) have been considered recently ||. Our results show that the 
model is not robust. Its main properties are radically changed, even for small disturbances 
in initial conditions and in the model itself, due to the appearance of real, gravitational 
singularities that are dynamically unavoidable in general. The singularities are, essentially, 
of two types. The first one corresponds to the hypersurfaces F(<p) = 0. It is not present 
in the isotropic case, and it implies that all previous homogeneous and isotropic solutions 
passing from the F((p) > to the F((p) < region are extremely unstable against anisotropic 
perturbations. The second type of singularity corresponds to F\{4>) — 0, with 

F 1 (<j>) = m + l(F'(<f>))\ (2) 

and it is present even for the homogeneous and isotropic cases. Although for small deviations 
of the conformal coupling the latter singularities are typically very far from the region of 
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interest, in the general case they can alter qualitatively the global dynamics of the model 
due to restrictions that it imposes on the phase space. Again, the persistence of some of our 
previously described results, in particular the ones concerning heteroclinic and homoclinic 
solutions, are challenged. 

Both kinds of singularities have already been described before. To the best of our knowl- 
edge, Starobinski f|J] was the first to identify the singularity corresponding to the hypersur- 
faces F(4>) = 0, for the case of conformally coupled anisotropic solutions. Futamase and 
co-workers || identified both singularities in the context of chaotic inflation in F(<f>) = 1— £0 2 
theories (See also The first singularity is always present for £ > and the second one 
for < £ < 1/6. Our conclusions are, however, more general since we treat the case of 
general F(<f>) and our results are based on the analysis of true geometrical invariants. Our 
main result is that the system governed by ([!]) is generically singular on both hypersurfaces 
F(4>) = and Fi((j>) = 0. Here, generically means that it is possible to construct non-singular 
models if one fine-tunes F(<p) and V(4>), as we will show below. 

One can advance that there are some geometrically special regions on the phase space 
of the model in question by a very simple analysis of the equations derived from the action 
(|I|). They are the Klein-Gordon equation 

+ V(0)it!=O, (3) 

and the Einstein equations 

F(<j>)G ab = (1 + F"(<f ) ))d a( f ) d b( f ) 

- l -g ab [(1 + 2F"(0))«9 C 0«9 C + 2V(<P)] - F'{4>) ( to D0 - V a 0V b 0) . (4) 

We will consider here the simplest anisotropic homogeneous cosmological model, the Bianchi 
type I, whose spatially flat metric is given by 

ds 2 = -dt 2 + a 2 {t)dx 2 + b 2 {t)dy 2 + c 2 (t)dz 2 . (5) 

The dynamically relevant quantities here are 
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Hi = -, H 2 = \, and H 3 = - . (6) 

a b c 

For such a metric and a homogeneous scalar field <fi = 4>(t), after using the Klein-Gordon 
Eq. (0), Eq. (f|) can be written as 

F(<f))G 00 = + V(<f>) - F'(<f)) {H, + H 2 + H 3 ) <P, (7) 
±F(d>)G u = 1 + 2 f W 2 - V(0 - Uj> + V'tf) - , (8) 



\ 2 F{4>)G 22 = 1 + 2 f /(0) 2 - V(0) - (tf 2 + V'(<j>) - , (9) 

- 2 mG^ = 1 + T (0 ^ 2 - - ^'(0) + V\4>) - ^r) . (10) 



c 2 



It is quite simple to show that Eqs. (|8|)-(|T0|) are not compatible, in general, on the hy- 
persurface F((f>) = 0. Subtracting @ and (|T0|) from ® we have, on such hypersurface, 
respectively, 

F'{4>){H\ - H 2 )<P = 0, and (H x - H 3 )<P = 0. (11) 

Hence, they cannot be fulfilled in general for anisotropic metrics. As it will be shown, it 
indeed corresponds to an unmovable (in the Painleve sense |7j]) geometrical singularity which 
cannot be prevented in general by requiring that F'((j)) = or <fi = on the hypersurface. 

As to the second singularity we have, after taking the trace of the Einstein equations, 
that: 

R = R{<f>, 0) = (4V(0) + W (0)F'(0) - (1 + F'(0))0 2 ) . (12) 

Inserting Eq. (|T2]) in the Klein- Gordon Eq. (^), one can see that it contains terms which are 
singular for Fi((p) = 0. Again, as we will see, this corresponds to an unmovable geometrical 
singularity, and it cannot be eliminated, in general, by demanding that F'((p) = on the 
hypersurface Fi(<f)) = 0. In both the hypersurfaces F(<p) = and Fi(<ft) = the Cauchy 
problem is ill-posed, since one cannot choose general initial conditions. 

The hypersurfaces F(<p) = and F\(<f>) = also prevent the global definition of an 
Einstein frame for the action (|l|), defined by the transformations 



g ab = F(<j))g ab , (13) 

(dJ\ 2 = FM_ 

It is well known that in the Einstein frame the Cauchy problem is well posed. Again, 
the impossibility of defining a global Einstein frame shed some doubts about the general 
Cauchy problem. Moreover, the standard perturbation theory for helicity-2 and helicity-0 
excitations, derived directly from Eqs. flT3|)-(|T4|), fails in both hypersurfaces ||. 

II. THE SINGULARITIES 

In order to check the geometrical nature of these singular hypersurfaces, let us consider 
the Einstein Eqs. ([7D-(|K]) i n detail. For the metric (||), we have the following identities 

Coo = H X H 2 + H 2 H 3 + H X H 3 , 

G u = a 2 (h, + H^H, + H 2 + H 3 ) - ±f?) , 

G22 = b 2 (H 2 + H 2 (Hx + H 2 + # 3 ) - ifl) , (15) 
G33 = c 2 (h 3 + H 3 (Hx + H 2 + H 3 ) - ~R 



R = 2[H 1 +H 2 + H 3 + H{ + H l 2 + H l 3 + H t H 2 + H 2 H 3 + H X H, 

After using expressions ( |15|) and introducing the new dynamical variables p = H\ + H 2 + H 3 , 
q = Hi — H 2 , and r = H\ — H 3 , Einstein Eqs. can be cast in the form 

1 d> 2 
E(cf>, 0, p, q, r) = --F{cj>) (p 2 + qr - q 2 - r 2 ) + ^ + V{(j>) - pF'{<j>)i> = 0, (16) 

r = -(p+^)r, (18) 

- 2FMP = (F(<t>) + 2F'(<P) 2 )p 2 + 2(1 + 2F"(0))0 2 - W{<j>) - 3F'(<f>)V'(cf>) 

-p<j>F'(4>) + (F((f>) + F'{(j)) 2 ){q 2 +p 2 - qr) (19) 

Using the energy constraint fll6"D , this last equation can be put in the form 



The Klein- Gordon equation (|3]) reads simply 

4> + p<P + V'(cf>) - ^-R((f>, <j>) = 0, (21) 

with -R(0, 0) given by Eq. (|T2|). The energy constraint, equation (|16D, is evidently compatible 
with the other ones. Indeed E((p, <p,p, q, r) — is an invariant surface since one has that 

j t E{<t>, fa p, q,r) = - (2p + ^^Aj E{4>, j>, P, q, r) (22) 

along solutions of Eqs. (0), (0), and fl2T|). Note that Eqs. (JHJ) and (0) are 



decoupled from the equations for q and r. Equations (Ti7| ) and ( jl£| ) are, hence, linear first 
order equations, and they could be easily integrated after the solutions of (|2H| ) and (^T|) have 
been found. Moreover, since one has rq — qr = 0, q(t)/r(t) is a constant of motion fixed 
only by the initial conditions. Suppose the initial ratio is g(0)/r(0) = 7: this would imply 
that (Hi — H 2 ) = J (Hi — H 3 ) for all t, leading to, for instance, c 7 (t) oc a'~ l (t)b(t) in the 
metric @. This simplification is a consequence of the scalar character of our source field, 
and it does not suppose any loss of generality in our arguments. 

A closer analysis of Eqs. (|17|)-(|2~0|) reveals the presence of the singularities. In general, 
the right-hand side of these equations diverge for F(<p) = and for Fi(<j)) = 0. One can check 
that these divergences are indeed related to real geometrical singularities by considering the 
Kretschman invariant / = R a t )C dR abcd , which for the metric @ is given by 

1 = 4 + Hi) 2 + (H 2 + Hi) 2 + (#3 + Hi) 2 + H\H\ + H\H\ + H 2 2 H^ . (23) 

As one can see, / is the sum of non negative terms. Moreover, any divergence of the variables 
Hi, H 2 , H 3 , or of their time derivatives, would suppose a divergence in /, characterizing a 
real geometrical singularity. Since the relation between the variables p, q, r, and Hi, H 2 , H$ 
is linear, any divergence of the first, or of their time derivative, will suppose a divergence in 
I. 

6 



Suppose, now, that F(<f>o) = 0, and that F((f>) is (real) analytical for = 0o. in this 
case, F'((f>)/F((f>) diverges as (0 — 0o) _1 near O for nonvanishing q and r, rendering q and 
r divergent by ( |T7]) and ([18]), unless p also diverges in order to keep (p + y^4>) finite on 
0o- In both cases, / diverges. There is no dynamical restriction to ensure that vanishes 
on O - it can take any value compatible with the energy constraint (|T6|). Indeed, the latter 
implies that on O 

^-pF'(0 o )0 + V(0 o ) = O. (24) 

There is no way of having = on O , unless V(4>o) = 0, and even in this case, = pF'((f> ) 
is also possible. Note that the hypothesis of F((f>) analytical at 0o is not a necessary one. For 
any different iable function F((f>) with a zero in O one has |-F(0)| = | jf F'(s)ds\ < k\<f) — <p \, 
with k = max se[ ^^] \F'(s)\, implying that \F'(<j))/F(<j))\ > \<f> - (po^F'^/k. Since F'(<j>) 
is assumed to be continuous, the last ratio tends to 1 when — > 0o, implying the divergence 
of F' ((f)) /F ((f)) in that limit. 

Now, let us suppose -Fi(0i) = 0. If F'((f>i) ^ 0, the right-hand side of Eq. ( pop diverges. 
The vanishing of F'((f>i) implies, by (^]), that F(<pi) = 0, and the arguments of the last 
paragraph can be repeated. 

A singularity-free model can be constructed by demanding that F((f) ) = F'((f> ) = 0, 
by choosing a V((f>) that goes to at a proper rate when — > (f> , and by demanding that 
F\((f>) have no other zeros than the ones of F((f>). Models for which F((f>) = £0 2n and 
V((p) = a<f) 2<y2n ~ lS} + high order terms, for instance, fulfill these requirements. However, such 
a highly fine-tuned class of model is of no physical interest here, since it does not contain 
F((f>) > and F((f>) < regions and consequently has no solution for which the effective 
gravitational constant G e R changes it sign along the cosmological history. The stability of 
such solutions were the starting point of the analyses of the pioneering work 0] and of the 
present one as well. Note that by Eq. (fj), models with an F((p) < region will allways 
have singularities of the type -F\(0) = 0. This fact shall be taken into account to better 
understand the recently proposed dynamical stability of the F((p) < region [||. 
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III. CONCLUSION 



The singularities described in the precedent section imply that the model presented in 
is not robust, since our main conclusions were a consequence of very especial initial 
conditions. For instance, all homogeneous and isotropic solutions crossing the F(4>) = 
hypersurface are extremely unstable against anisotropic perturbations. By Eqs. flT7| ) and 
(Pf), any deviation from perfect isotropy (expressed by nonvanishing q and r variables) for 
these solutions, however small, will lead catastrophically to a geometrical singularity. Many 
of the novel dynamical behaviors presented in |TJ,|2| depend on these solutions. This is the 
case, for instance, of some solutions exhibiting superinflation regimes. The heteroclinic and 
homoclinic solutions identified in |],^|] can cross the F(<f)) = hypersurface and, hence, they 
also suffer the same instability against anisotropic perturbations. The homoclinic solutions 
were considered as candidates to describe a non-singular cosmological history, with the 
big-bang singularity being avoided through a classical birth of the universe from empty 
Minkowski space. Apart from F((j>) = singularities, these solutions are also affected by 
the singularities of type F\{4>) — 0- Suppose that the conformal coupling is disturbed by a 
very small negative term: F(<p) = 1 — (| — e)</> 2 . The F 1 (<p) = singularities will be near 
the <p = H/y/e hypersurfaces. Although they are located far from the F(<f)) = regions, 
they alter the global structure of the phase-space. In this case, they restrict the existence 
of homoclinics, rendering a non-singular cosmological history more improbable. 

The singularities do not affect the conclusions obtained by considering solutions inside 
the F(<p) > region. The asymptotic solutions presented in f2|, for instance, are still valid. 
The conclusion that for large t the dynamics of any solution (inside F((p) > 0) tends to an 
infinite diluted matter dominated universe remains valid. Moreover, for small anisotropic 
deviations (q and r small in comparision with p), Eqs. (|17|) and (|18|) allow us to conclude that 
solutions inside F((p) > 0, for large t, approach exponentially isotropic matter-dominated 
universe. 
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